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gf  SUMMARY 

A three-dimensional  steady  state  mathematical  model  of  a cable  used  to  suspend 
a sonar  body  from  a helicopter  is  derived.  Results  are  presented  in  dimensionless  form  for 
a length  of  cable  totally  immersed  in  any  Newtonian  fluid,  and  may  therefore  be  applied 
to  a number  of  practical  problems.  Various  approximations  are  made  in  order  to  obtain 
analytical  solutions.  The  suitability  of  the  approximations  is  investigated  by  comparing 
solutions  of  the  general  and  approximate  equations  for  a uniform  fluid  velocity.  In  general, 
where  maximum  angular  deviations  from  the  vertical  are  small,  agreement  is  observed 
to  be  very  good.  Two-dimensional  i power * and  llogl  laws  are  assumed  for  the  wind 
above  the  sea,  and  a three-dimensional  Ekman  spiral  is  assumed  for  the  sea  current. 
For  these  non-uniform  fluid  velocities,  the  concept  of  an  ♦ equivalent  mean  fluid  velocity ♦ 
is  used  to  simplify  greatly  the  presentation  of  results. 


POSTAL  ADDRESS:  Chief  Superintendent,  Aeronautical  Research  Laboratories, 
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1.  INTRODUCTION 
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For  the  underwater  detection  of  submarines,  helicopters  commonly  use  a sonar  body  which 
is  lowered  by  cable  into  the  water.  For  efficient  operation  of  the  system,  it  is  generally  important 
for  the  sonar  body  to  remain  as  close  to  the  vertical  as  possible  while  submerged.  For  a calm 
sea  with  little  or  no  wind,  this  is  achieved  by  manoeuvring  the  helicopter,  either  manually  or 
automatically  via  the  helicopter  control  systems,  in  such  a way  that  the  cable  deviations  (in 
pitch  and  roll)  from  the  vertical  at  the  suspension  point  are  kept  close  to  zero.  For  the  case  of  a 
steady  non-zero  wind,  resulting  cable  bow  requires  one  or  both  of  these  angular  deviations  to  be 
kept  constant  at  non-zero  values.  A mathematical  model  of  the  cable  and  sonar  body  enables 
these  angles,  which  are  referred  to  as  “offset  cable  angle  corrections”,  to  be  estimated  and  allows 
an  investigation  of  the  behaviour  of  the  system  under  certain  conditions.  While  hovering  above 
the  sea,  the  helicopter  normally  heads  into  wind  and  both  speed  of  helicopter  relative  to  ground 
and  wind  speed  relative  to  helicopter  are  usually  given  by  the  helicopter  instrumentation.  It  is 
assumed  therefore  that  these  quantities  are  given  and  that  no  further  data  about  the  air  and  sea 
environment  are  measured  or  known.  A mathematical  representation  of  the  air  and  sea  must 
therefore  be  assumed  in  order  to  calculate  the  fluid  drag  force  on  the  cable  and  sonar  body. 

A three-dimensional  time-dependent  model  in  which  the  cable  is  represented  by  a number 
of  rigid  links  is  described  in  Reference  I . However,  it  is  sufficient,  and  considerably  more  con- 
venient, to  consider  a steady  state  model  to  investigate  much  of  the  behaviour  of  the  system. 
This  applies  particularly  to  the  problem  of  calculating  the  offset  cable  angle  corrections,  which 
need  to  be  assumed  invariant  with  time.  The  treatment  ignores  the  sea  waves  which  would  cause 
non-steady  effects,  including  cable  oscillation.  Steady  state  solutions  have  been  presented  by 
other  authors,-'3'1  but  are  restricted  to  two-dimensional  motion  in  a uniform  flowfield.  These 
restrictions  are  removed  in  this  document.  Results  are  presented  in  dimensionless  form  for  a 
length  of  cable  totally  immersed  in  any  Newtonian  fluid,  and  may  therefore  be  applied  to  a 
number  of  cable  problems. 

The  general  equilibrium  equations  are  first  derived  in  Section  2.  Solutions  to  these  equations 
giving  cable  orientation  at  the  top  end  of  the  cable  are  then  presented  in  Section  3.  The  equations 
are  solved  numerically  using  the  Runge-Kutta  second  order  method  for  the  general  case  and 
analytically  for  the  two-dimensional  case  with  uniform  fluid  velocity. t Various  approximations 
are  then  made  which  allow  an  analytical  solution  for  the  three-dimensional  case.  The  suitability 
of  the  approximations  is  shown  by  parameter  studies  in  which  uniform  fluid  velocity  solutions 
of  the  general  and  approximate  equations  are  compared. 

In  Section  4,  it  is  shown  that  a non-uniform  fluid  velocity  may  be  replaced  by  an  “equivalent 
mean  fluid  velocity”,  thus  greatly  simplifying  the  presentation  of  results.  This  concept  is  also  of 
use  in  developing  a mathematical  model  for  a flight  simulator,  where  it  is  important  to  simplify 
the  representation  of  the  environment  as  much  as  possible.  Equivalent  mean  fluid  velocities  are 
then  presented  for  typical  v/ind  and  sea  velocity  profiles.  These  include  two-dimensional  “power" 
and  “log”  laws  for  the  wind  above  the  sea,  and  a three-dimensional  Ekman  spiral  for  the  sea 
current.  Comparisons,  similar  to  those  for  the  uniform  fluid  velocity  case,  are  made  between 
solutions  of  the  general  and  approximate  equations.  Finally,  an  example  is  given  that  is  typical 
of  a sonar  body  suspended  from  a helicopter. 

2.  MATHEMATICAL  MODEL 
2.1  Approach 

In  order  to  investigate  the  behaviour  of  a cable  used  to  suspend  a sonar  body  from  a heli- 
copter, the  following  steady  state  mathematical  model  is  formulated.  Given  the  tension  and 

t Derivations  of  analytical  solutions  are  given  in  Reference  2 but  are  repeated  here  for 
completeness  of  presentation  in  the  notation  of  this  document. 


» 
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angular  orientation  of  the  cable  at  the  attached  body,  the  angular  orientation  of  the  cable  at  the 
helicopter  is  required  for  any  given  steady  state  relative  motion  between  the  helicopter  and  the 
air  and  sea.  When  the  body  is  submerged,  a solution  is  first  obtained  for  the  length  of  cable  also 
submerged.  This  solution  then  provides  the  necessary  boundary  conditions  for  the  solution  for 
the  length  of  cable  above  the  sea. 

Consider  then  a length  of  cable,  totally  immersed  in  a Newtonian  fluid,  with  the  fluid  moving 
horizontally.  It  is  assumed  that  the  only  forces  acting  on  the  cable  are  the  tension  force,  effective 
gravitational  force  (i.e.  including  buoyancy  effect),  and  fluid  drag  force  normal  to  the  cable. 

Effects  of  cable  twist,  skin  friction,  induced  mass,  and  elasticity  of  the  cable  are  not  considered 
here.  The  tension  force  and  angular  orientation  at  P i,  the  lower  end  of  the  cable  length,  are 
assumed  given.  For  steady  state  motion,  the  angular  orientation  at  Po,  the  upper  end  of  the  cable 
length,  is  required. 

2.2  Axes  Systems 

It  is  convenient  to  define  the  following  four  right-handed  orthogonal  axes  systems: 

(a)  Cable  axes — xyz 

The  origin  is  at  Po-  The  z axis  is  vertical  and  positive  downwards.  The  x axis  is 
horizontal  and  lies  in  the  fuselage  longitudinal  plane  of  symmetry  (positive  forward). 

The  y axis  forms  a right-handed  orthogonal  axes  system. 

(b)  Earth  axes—xey’BZE 

The  earth  axes  coincide  with  the  cable  axes  at  any  arbitrarily  chosen  instant  of  time. 

(c)  Cable  element  axes — x,y,z, 

Cable  element  axes  are  defined  for  each  element  of  length  8s  of  cable.  The  origin  P 
is  at  a distance  s along  the  cable  from  Po  and  has  coordinates  (x,  y,  z)  in  cable  axes. 

The  zs  axis  is  tangential  to  the  cable  at  P and  is  positive  in  the  normally  downward 
direction.  The  x,  axis  is  normal  to  the  z,  axis  and  lies  in  the  vertical  plane  containing 
the  z,  axis.  The  y,  axis  forms  a right-handed  orthogonal  axes  system. 

(d)  Intermediate  rotation  axes — x'y'z' 

Intermediate  axes  are  obtained  when  rotating  xyz  axes  to  give  x»y,z,  axes  (see  below 
and  refer  to  Figure  1). 

Earth  axes  serve  only  to  define  the  separate  cable  and  fluid  motions  with  respect  to  earth. 

For  the  motion  of  fluid  relative  to  the  cable  required  in  the  analysis  below,  cable  axes  are  used. 

It  should  be  noted  though  that  because  of  the  invariance  with  time  of  the  cable  and  fluid  motion, 
the  xf. , ys,  and  ze  directions  are  always  the  same  as  the  respective  x,  y,  and  z directions. 

Angular  orientation  of  the  cable  element  axes  relative  to  the  cable  axes  is  expressed  in  terms 
of  the  conventional  Euler  angles  used  in  aircraft  aerodynamics  where  the  order  of  rotation  is  1 

yaw,  pitch,  and  roll.  However,  because  effects  of  cable  twist  are  ignored,  only  two  angles  are 
required  to  define  the  orientation.  In  general,  for  the  conditions  assumed,  yaw  angle  will  vary 
between  — 1 80°  and  1 80°,  whereas  the  magnitude  of  pitch  and  roll  angle  will  always  be  less  than 
90°.  Hence,  by  choosing  to  use  pitch  and  roll  angle,  singularities  that  occur  at  angles  of  ±90° 
are  avoided.  Following  parallel  translation  of  the  cable  axes  such  that  their  origin  coincides 
with  the  origin  of  the  cable  element  axes,  Figure  1 shows  the  appropriate  angular  rotations. 

The  xyz  axes  are  first  rotated  through  an  angle  8 about  the  y axis  to  give  intermediate  axes, 
x'y’z'.  These  axes  are  then  rotated  through  an  angle  <f>  about  the  x'  axis  to  give  the  cable  element 
( axes,  x,y,z,. 

Although  the  pitch  and  roll  angles,  8 and  <f>,  are  mostly  used  in  the  following  analysis,  in 
some  instances  it  is  convenient  to  use  angles  obtained  by  a rotation  in  yaw,  y,  followed  by  pitch, 
y.  The  angles  may  be  related  to  8 and  <f>  by  considering  the  direction  cosines  of  the  cable  element 
axes  with  respect  to  the  cable  axes.  The  direction  cosines  are  given  by  Duncan5  as 

cos  « = dx/ds  = sin  8 cos  <f> 

• cos  /?  = dy/ds  = —sin  <f> 

cos  y = dz/ds  = cos  8 cos  <f> 

where  the  pitch  angle  y is  given  directly  by  Equation  (lc),  and  the  yaw  angle  y is  given  by 

tan  y = cos  /3/cos  a = —cosec  8 tan  <f> 


(la) 
(1  b) 
(lc) 

(2) 
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2.3  Matrix  Transformations  for  Rotations  Through  Pitch  and  Roll  Angles 


To  resolve  the  various  forces  into  the  one  axes  system,  it  is  convenient  to  first  present  some 
useful  matrix  transformations.6 

Rotation  of  a vector  through  angles  0 and  <!>  is  performed  by  premultiplying  the  vector  by 
matrices  0 and  <l>  respectively,  where 


cos  0 

0 

— sin  0' 

1 

0 

0 ‘ 

0 

0 

1 

0 

, <I>  = 

0 

cos  <j> 

sin  </> 

sin  0 

0 

cos  0 

0 

—sin  <f> 

cos 

and  the  respective  matrix  inverses  are 


cos  0 

0 

sin  0' 
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0 ' 
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0 
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, <D  1 = 
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COS  i f> 

—sin  <f> 

—sin  0 
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cos  0 
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The  matrix  transformation  obtained  by  rotating  through  an  angle  0 followed  by  a rotation 
through  <t>  is 


<D0 


cos  0 0 

sin  0 sin  cos 

sin  0 cos  <j>  — sin  ^ 


— sin  0 
cos  6 sin  <j> 
cos  0 cos  <j> 


with  (he  inverse  given  by 


cos  0 

sin  0 sin  </> 

sin  0 cos  <j> 

[<I>  (-)]  i = 0 i <J>  i 

0 

COS  c/> 

—sin  $ 

—sin  0 

cos  0 sin  <j> 

cos  0 COS  i f> 

(5) 


(6) 


Inspection  shows  that  the  rotation  matrices  have  the  property  that  their  inverses  are  the 
same  as  their  transposes. 


2.4  Velocity  of  Cable  Relative  to  Fluid 

It  is  assumed  that  both  the  cable  and  fluid  may  move  in  any  horizontal  direction  defined 
by  the  angle  between  the  velocity  direction  and  the  x direction,  the  angle  increasing  positively 
in  the  usual  right-handed  sense.  The  magnitude  of  the  cable  and  fluid  velocities  relative  to  earth 
are  defined  as  uh*  and  uy*  respectively,  with  corresponding  directions  <Ph  and  i/>f.  The  quantities 
uh*  and  i pH  are  assumed  uniform,  and  uy*  and  >py  are  assumed  to  be  known  functions  of  z. 
The  magnitude  and  direction  of  the  velocity  of  the  fluid  relative  to  the  cable  are  defined  as  u* 
and  </>,  with  x,  y,  and  z components  given  by  the  vector 


u 

Uy 

uh 

V 

= 

Vy 

— 

vh 

w 

0 

0 

where 


(7) 


14 y = Uy*  COS  ipy, 

II h = Un*  cos  4>h, 
u = u*  cos  i fi. 


Vy  = Uy*  sin  ipy 
VH  = Uh*  sin  <pH 
v = u*  sin  </> 


On  resolving  in  cable  element  axes,  the  relative  fluid  velocity  vector  is  given  by 


u,' 

V 

U COS  0 

V, 

<1)0 

V 

= 

u sin  0 sin  <f>  + v cos  ^ 

w. 

0 

u sin  0 cos  <f>  — v sin  <f> 

The  velocity  component  w,  is  not  required  since  skin  friction  is  ignored.  The  velocity  magnitude 
and  direction,  u„*  and  </>,,  in  cable  clement  axes,  are  defined  by 


u,*  = (u,2  + v,2)* 
<P,  = arctan  (v,/u,) 


(9a) 

(9b) 


2.5  Equilibrium  Equations 

For  a cable  element  of  length  8s  at  P (see  Figure  2),  the  gravitational,  fluid,  and  tension 
forces  acting  on  the  element  are  first  resolved  in  cable  element  axes.  Since  the  relative  motion  is 
steady  state,  these  forces  are  in  equilibrium  so  that 

Tension  Force  + Fluid  Force  + Gravitational  Force  = 0 
which  may  be  rewritten  using  corresponding  vectors  expressed  in  cable  element  axes  as 


(10) 


The  gravitational  and  fluid  forces  are  assumed  to  act  at  P,  the  upper  end  of  the  cable 
element.  The  gravitational  force  vector  is  given  by 


Tx' 

Ty 

+ 

Fx 

Fy 

+ 

Gx 

Gy 

O' 

0 

Tz 

Fz. 

Gz 

0 

Gx 

' 0 ' 

—sin  0 

Gy 

= ct>0 

0 

= mg  8s 

cos  0 sin  <f> 

Gz 

mg  8s 

cos  0 cos  <f> 

(11) 


where  m is  the  effective  cable  mass  per  unit  length  (allowing  for  buoyancy  effect)  and  g is  the 
gravitational  acceleration  (along  the  z axis).  The  fluid  force  vector  is  given  by 


(12) 


Fx 

'cos  ip,' 

Fy 

Fz 

= kc  8s  (a,*)2 

sin  ip, 
0 

= kc  8s  u* 

s?  O 

i 

where 


kc  = IpF  dc  Cd 


(13) 


and  pf  is  the  fluid  density,  dc  is  the  cable  diameter,  and  Cd  is  the  drag  coefficient  of  the  cable  in 
the  lateral  direction  of  flow. 

Let  the  tension  force  at  P acting  on  the  upper  end  of  the  element  be  T (i.e.  along  z,  axis  in 
negative  direction).  Then  the  tension  force  at  the  lower  end  is  T + 87  with  orientation  relative 
to  cable  axes  given  by  a rotation  in  pitch  of  0 + 80  followed  by  a rotation  in  roll  of  <f>  + 8<£.  The 
resultant  tension  force  vector,  in  cable  element  axes,  acting  on  the  element  is  then  given  by 


(14) 


where  the  matrices  0t  and  <b,  are  used  to  rotate  a vector  through  angles  0 + 86  and  <£  + 8 <p 
respectively,  and  are  defined  by  substituting  0 + 80  for  0,  and  <f>  + 8<p  for  <p,  in  Equations  (3). 
Similar  substitution  in  Equation  (6)  defines  ©r1  ®i“l.  On  neglecting  second  order  and  higher 
order  differences  so  that 

sin  (0  + 80)  = sin  0 + 80  cos  0 (15a) 

cos  (0  + 80)  = cos  0 — 80  sin  0 (156) 

with  similar  relations  for  sin  (<p  + 8<p)  and  cos  (<p  + 8<p),  Equation  (14)  may  be  simplified  to 
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1 

- 7 

0 

1 

— 

- 7 8<p 
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On  substituting  Equations  (11),  (12),  and  (16)  in  Equation  (10), 


'7  80  cos  i f 

V 

— sin  0 

O' 

~T8<p 

+ kc  8s  u,* 

v. 

+ mg  8s 

cos  0 sin  < p 

= 

0 

8T 

.0. 

cos  0 cos  <p 

0 

(16) 


(17) 


In  the  limit,  as  8s  tends  to  zero,  Equation  (17)  may  be  written  in  differential  form  as 

4 


(18a) 


w 


= mg  sin  8 — kc  u,*  u , 

= mg  cos  6 sin  <f>  + kc  u,+  v,  (18ft) 

— mg  cos  9 cos  <(>  (18c) 

as 

Using  Equation  (lc),  where  dz/ds  = cos  9 cos  <f>,  the  above  equations  may  be  transformed  to 


dO  __ 

mg  sin  8 — kc  us*  u, 

(19a) 

dz 

T cos  8 cos2  4> 

d<f> 

mg  cos  8 sin  <f>  + kc  u,*  v. 

T cos  8 cos  <l> 

U **>) 

dT 

dz 

- mg 

(19  c) 

which  are  the  basic  equilibrium  equations,  in  differential  form,  used  here. 


T cos  <f> 


dO 

ds 


ds 


3.  SOLUTION  OF  EQUILIPRIUM  EQUATIONS 

In  their  most  general  form,  Equations  (19a,  ft,  c)can  only  be  solved  numerically.  However, 
analytical  solutions  are  possible  given  further  assumptions  or  approximations.  Solutions  of  both 
the  general  and  approximate  equations,  which  are  termed  “general”  and  “approximate”  solutions 
respectively,  are  given  below.  As  a convenient  reference,  a classification  of  these  solutions  is 
given  in  Table  1. 


TABLE  1 

Classification  of  Types  of  Solution  of  Equilibrium  Equations 


Solution  type 

Method  of  solution 

1.  General  (i)  2-  or  3-D  variable  fluid 
velocity 

Numerical  solution  of  Equations  (22a,  b) 
using  Runge-Kutta  2nd  order  method 

(ii)  Restricted  to  2-D  uniform 
fluid  velocity 

Analytical  solution  given  by  Equation  (26) 

2.  Approx.*  (i)  2-  or  3-D  variable  fluid 
velocity 

Approximate  solution  given  by  Equations 
(33a,  ft)  requiring  evaluations  of  integral 
by  Simpson’s  formula 

(ii)  Restricted  to  uniform  fluid 
velocity! 

Analytical  solution  given  by  Equations 
(34a,  ft) 

3.1  General  Equations 

To  obtain  a solution  of  Equations  (19a,  ft,  c),  Equation  (19c)  is  first  integrated  between  P 
and  Pi  (where  z = z\  and  T = T\),  to  give 


» T , 

dT  = 
J T 


(20) 


* Approximations  are  (a)  cable  mass  concentrated  at  attached  body  and  (ft)  small  angle 
approximations. 

t Including  non-uniform  fluid  velocity  if  equivalent  mean  fluid  velocity  components  are 
obtained  analytically. 


' 
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which  results  in 


i 


> 


T = Ti  - mg  (z  - z i)  (21) 

The  dimensionless  quantities  /x  = mgzi/Ti,  X = x/zi,  T = yi/zi,  Z = z/zi,  5 = s/zi,  and 
t/*  = (Arc  zi/Ti)*u*  (similarly  for  all  velocities,  e.g.  V = (ke  zi/7i)*v)  are  now  introduced  with 
the  subscripts  already  used  retaining  their  meaning.  On  substituting  for  T from  Equation  (21), 
Equations  (19a,  b)  may  be  written  in  terms  of  the  dimensionless  quantities  as 

d8  _ p sin  8 - u.*  u, 

dZ~  [ 1 ~/x(Z  - 1)]  cos  flcos2  4>  (22a) 

d<\>  /x  cos  8 sin  + {/,*  K, 

</z  = [ i - #x(z  - l)]  cos  e cos  t (22A) 


Given  0 = 0i  and  <f>  = <£i  at  Pi,  the  above  pair  of  simultaneous  linear  differential  equations 
must  be  solved  in  order  to  obtain  the  angles  8 — 80  and  <j>  = <f>o  at  P0.  where  Z = 0.  Knowing 
the  angles  8 and  <f>,  the  cable  shape  can  readily  be  evaluated  by  numerical  integration  of  Equations 
(la,  b,  c)  (see  Part  1 of  Appendix  A). 

3.1.1  Numerical  Solution  for  General  Case 


For  the  general  case,  Equations  (22a,  b)  are  solved  numerically  using  a Runge-Kutta  second 
order  method.  For  the  special  case  described  in  Section  3.1.2,  i.e.  two-dimensional  case  with  a 
uniform  fluid  velocity,  an  analytical  solution  is  possible.  A direct  check  on  the  accuracy  of  the 

numerical  method  may  therefore  be  made  for  this  case.  This  was  done  to  first  check  the  method  \ 

and  then  to  obtain  a suitable  integration  interval  to  be  used  for  more  general  conditions  where 

an  analytical  solution  is  not  possible.  When  used  for  each  different  type  of  condition,  further 

checks  were  made  in  the  usual  way  by  comparing  solutions  for  varying  integration  intervals.  It 

was  found  that  an  integration  interval  in  Z of  0- 02  gave  sufficient  plotting  accuracy  for  the 

curves  presented  in  the  figures. 


3.1.2  Analytical  Solution  for  Special  Case 

For  the  special  case  of  a two-dimensional  solution  for  a uniform  fluid  velocity,  an  analytical 
solution  of  Equations  (22a,  b)  is  possible. t A two-dimensional  solution  in  any  arbitrarily  rotated 
vertical  plane  is  most  directly  given  by  using,  instead  of  8 and  <f>,  the  variables  x and  y,  where 
Ixl  is  constant  and  y > 0.  However,  because  of  possible  sign  changes  in  x,  it  is  more  convenient 
to  consider  the  xz  plane,  using  8 and  <f>,  where  8 may  be  positive  or  negative,  and  <f>  = 0.  For  a 
subsequent  rotation  of  the  vertical  plane,  y is  then  given  by  the  value  of  \8\  previously  calculated 
for  the  xz  plane. 

Setting  <l>  = 0,  V = 0,  and  U = Uo,  then  U,*  = |C/»|,  U,  = Uo  cos  8,  and  Equations  (22a,  b ) 
reduce  to  the  single  equation 

d8  /x  sin  8 — Uo  | C/o|  cos2  8 

dZ  = [1  - fx(Z  — 1)]  cos  8 ( ) 


On  integrating  Equation  (23)  between  Po  and  P, 


P * 

dZ 

U-KZ-Dl 

J z 

Putting  l = sin  8,  a = /x/[t/0  |l/o|],  q = (1 
be  written  as 


[fj.  sin  8 


+ 4/a2)*,  and 


cosO dg 

- Uo  |f/0|  cos2  8} 


(24) 


b = 1 + /x(l  — Z),  Equation  (24)  may 


t Refer  to  footnote  on  page  1. 


i 


6 


( 
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In  b = a 


dt 

(»*  + at  - 1) 


[2/  - a(q  - 1)]| 

[2/  + a(q  + 1)]J 

t 

__  1 |[2/i  — a(q  — 1)]  [2 / + a(q  + l)]) 

q |[2/i  -f  a(<7  + l)J  [2/  ~ a(q  - 1)]J 

which  gives 


. 2(6«  - 1)  + a[b*{q  - 1)  + (?  + 1)]  sin  0i 

sin  v = — — — — 

a[(q  - 1)  + b«(q  + 1)]  + 2 (bo  - 1)  sin 

Asa  special  case,  when  Si  = 0,  Equation  (26)  becomes 


sin  8 = 


2 (b<>  - 1) 

a[(q  - 1)  + b«(.q  + 1)] 


The  angle  8 = 80  at  Po  is  given  by  setting  Z = 0 (i.e.  b = l + n). 


(25) 


(26) 


(27) 


3.2  Approximate  Equations 

To  obtain  an  analytical  three-dimensional  solution,  the  following  approximations  are  found 
to  be  convenient : 

(a)  the  “effective”  cable  mass  is  considered  to  be  concentrated  at  the  attached  body; 

(b)  small  angle  approximations  for  8 and  are  assumed  (i.e.  sin  8 x 8,  cos  8 x 1,  and 
similarly  for  <£). 

For  approximation  (a)  to  be  valid,  either  the  ratio  of  cable  mass  to  body  mass,  or  the 
angles  8 and  <f>,  must  be  small. 

Approximation  (a)  is  achieved  by  setting  m = 0 in  Equations  (19a,  b,  c)  to  give 


dd  — kc  us*  u, 

(28a) 

dz  T cos  8 cos2  <f> 

d<f>  kc  u,*  v. 

(286) 

dz  T cos  8 cos  <!> 

dT  „ 

<fe=° 

(28c) 

where  the  tension  T\  at  P\  is  now  replaced  by  7\  + mgzi.  From  Equation  (28c),  the  tension 
along  the  cable  is  seen  to  be  constant  and  is  given  by 

T = fi  + mgzi  = 7-1(1  + m)  (29) 

Substituting  for  T from  Equation  (29),  Equations  (28a,  b ) may  be  written  in  dimensionless 
form  as 


dd  - U.*  U. 

dZ,  (1  + ft)  cos  8 cos2  <f> 

d<f>  U,*  V, 

dZ  (1  + fj.)  cos  8 cos  4> 


(30a) 

(306) 


Using  approximation  (6),  then  from  Equation  (8),  U,  x U,  V,  x V,  and  U ,*  x U*,  and 
Equations  (30a,  b ) may  be  written  as 


dd  - U*  U 
dZ  (1  +h) 

d<f,  U*  V 
dZ  (1  + h) 


(31a) 

(316) 
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In  Equations  (31a,  b)  and  below,  the  variables  U and  V may  be  alternatively  expressed  in  terms 
of  U*  and  <fi  using  the  relations 


(/=  U*  cos  ip 

(32 a) 

V=  U*  sin^- 

(326) 

On  integrating  Equations  (31a,  b)  between  P and  Pi, 

8 =8i  + 1 U*  UdZ 

(1  + P) 

(33a) 

J z 

/%  | 

*=* U*VdZ 

+ ^ J z 

(33  b) 

where  the  angles  8 80  and  <f>  <j> o at  P0  are  given  by  setting  Z = 0.  Equations  (33a,  b)  provide 

an  analytical  three-dimensional  solution  provided  a suitable  fluid  velocity  function  is  chosen 
such  that  the  integrations  can  be  performed  analytically.  The  cable  shape  may  readily  be  obtained 
using  the  above  small  angle  approximations  (see  Part  2 of  Appendix  A). 

3.3  Comparison  Between  Uniform  Fluid  Velocity  Solutions  of  General  and  Approximate 
Equations 

For  a uniform  fluid  velocity  (i.e.  U*  = Uo*,  U = Uo,  V = Vo,  and  <ji  = < !>0)  the  integrations 


in  Equations  (33a,  b)  can  be  performed  analytically  to  give 

8 = «i  + Uo * Uo  (1  - Z)/(l  + /i)  (34a) 

<f>  = <h  ~ Uo * Vo  (1  - Z)/(l  + p)  (346) 

or,  using  Equations  (32a,  b), 

8 = 8 , + (Uo*)2  cos  0O(1  - Z)/(l  + ti)  (3 Sa) 

(W)2  sin  0od  - Z)/(l  + p)  (356) 


The  validity  of  the  “approximate”  solutions  is  discussed  below  by  comparing  them  with 
appropriate  “general”  solutions  for  both  two-dimensional  and  three-dimensional  cases. 


3.3.1  Two-dimensional  Solution 

For  the  two-dimensional  case  with  uniform  fluid  velocity,  the  “general”  solution  is  given 
analytically  by  Equation  (26).  To  obtain  a corresponding  “approximate”  solution,  set  Uo*  = |Uo| 
in  Equation  (34a).  Hence,  the  “approximate”  solution  in  the  xz  plane  is 

8 = fli  + Uo  |Uo|(l  - Z)/(l  + n)  (36) 

At  Z — 0,  comparisons  between  this  solution  and  the  “general”  solution  are  given  in  Figures  3 
to  5 for  appropriate  values  of  the  parameters  p,  8\,  and  Uo,  where  Uo  > 0.  For  Uo  < 0,  symmetry 
considerations  readily  provide  appropriate  solutions,  i.e.  sign  changes  in  Uo  and  8 1 result  in  a 
sign  change  in  do.  When  considering  the  length  of  cable  in  air,  negative  Uo  is  the  more  usual 
case  since  the  helicopter  normally  heads  into  wind.  A few  results,  which  show  cable  shapes  for 
various  parameter  values,  are  presented  in  Figure  6.  It  can  be  seen  in  Figures  3 to  6 that  agree- 
ment is  good  where  maximum  angular  deviations  of  the  cable  from  the  vertical  are  small  (usually 
for  small  values  of  the  parameters  n,  8it  and  Uo).  As  one  would  expect  when  these  deviations 
become  large,  approximation  ( b ),  and  hence  approximation  (a)  when  n is  not  small,  become 
increasingly  less  valid,  resulting  in  poor  agreement. 

3.3.2  Three-dimensional  Solution 

In  presenting  three-dimensional  solutions,  it  is  convenient  to  specify  the  boundary  con- 
ditions at  Pi  in  terms  of  *i  and  yi,  instead  of  8 1 and  <f> i.  From  Equations  (lc)  and  (2),  8 and  ^ 
are  given  directly  from  the  relations 
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tan  0 = tan  y cos  \ (37fl) 

sin  <j>  = - sin  y sin  x (37 b) 

b so  that  Oi  and  <f>\  may  be  readily  obtained.  At  Z = 0,  Figures  7 to  9 show  comparisons  between 

three-dimensional  solutions  of  angular  orientation  in  pitch  and  roll  obtained  using  the 
“approximate”  solutions  given  by  Equations  (35a,  6)  and  the  “general”  solution  of  Equations 
(22a,  b)  for  appropriate  values  of  the  parameters  /*,  *i,  yi,  Uq*,  and  ifio.  For  three  specific  cases, 
corresponding  comparisons  between  cable  shapes  are  illustrated  in  Figure  10.  In  these  figures, 
where  an  arbitrary  choice  of  i/j0  is  required,  the  value  ±180°  is  chosen;  as  previously  stated,  this 
corresponds  to  the  usual  case  for  the  length  of  cable  in  air.  As  found  with  the  two-dimensional 
case,  where  maximum  angular  deviations  from  the  vertical  are  small  (usually  for  small  values  of 
the  parameters  /x,  yu  and  U0*),  agreement  is  observed  to  be  good;  large  deviations  result  in  poor 
agreement. 

4.  SOLUTIONS  FOR  TYPICAL  WIND  PROFILES  AND  SEA  CURRENTS 

For  the  case  of  a non-uniform  fluid  velocity,  it  is  convenient  to  obtain,  as  a first  step, 
“equivalent  values”  of  the  mean  fluid  velocity  components  Uf  and  VF.  They  are  defined  to  be 
those  values  which  give  the  same  solution  of  Equations  (33 a,  b)  at  Z = 0 as  obtained  when  the 
actual  non-uniform  values  are  used.  This  requires 


<t/*>  <C/>  = U * VdZ 

(38a) 

J 0 
/%  | 

<6/*><K>  = U*  VdZ 

(386) 

J 0 

where  all  variables  enclosed  by  “<”  and  “>”  are  mean  values.  For  any  non-uniform  velocity, 
once  these  equivalent  values  are  calculated,  solutions  giving  angular  orientation  may  be  readily 
obtained  using  the  uniform  velocity  results  of  Section  3.3.  This  procedure  greatly  simplifies  the 
presentation  or  results. 

At  this  stage,  non-dimensional  velocities  are  scaled  by  Ufo  * (assumed  non-zero),  the  value 
of  Uf*  at  Po , to  give  u>  — Uf*/Ufo*,  A Uu*/Ufo*,  and  a = U*/Ufo*  with  x and  y com- 
ponents defined  as 

<*>z  — w cos  ipF  Uf/Ufo  *,  OJy  — w sin  </>f  = V f/Ufo  * 

A*  = A cos  ipu  = Uh/Ufo  *,  Xy  = X sin  1/0/  = Vh/Ufo  * 

az  = a cos  ip  — U/Ufo  *,  </„  = a sin  ip  = V/Ufo  * 

Assuming  the  scaled  non-dimensional  fluid  velocity  components  wx  and  cu„  are  known  functions 
of  Z,  equivalent  mean  values  <aiz>  and  (wy)  are  now  derived.  Corresponding  non-scaled  values 
(Uf/  and  /Vf / may  be  obtained  by  multiplying  by  Ufo  *■  From  Equations  (7)  and  (32a,  b ),  the 
variables  in  Equations  (38a,  b)  may  be  expressed  in  terms  of  the  above  scaled  velocity  components 


as 

U — Uf  — Uh  = Ufo*  °z  (39a) 

V — VF  — Vh  — Ufo*  (39 b) 

U * = (U*  + F2)*  = | £/«,*|  (ax2  + ay2)*  (39c) 

<£7>  = Ufo*  <a*>  (40a) 

<F>  = Ufo*<*w>  (406) 

(U*y  = |LVo*|  «a*>2  + <ay>2)*  (40c) 

where  <oz/  - (a>z/  — Xz  and  (ay)  = (wy)  — Xy. 


Substituting  Equations  (39 a,  b,  c)  and  (40a,  b , c)  in  Equations  (38a,  6), 

«a*>2  + <a„>2)*  <a*>  = flz 
«a*>2  + <o„>2)*  <av>  = fty 
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(41a) 

(416) 


U*  UdZ 


(oz  2 + fv2)'  oz  dZ 


U*VdZ  = (oz2  + Oy2)*  ay  dZ 


Equations  (4 1 a,  b)  may  be  solved  for  /tuz')  and  <<o„>  to  give 


‘X?  A x + «z  \Px\/Px 


— Aj,  4-  <*y 


for  f}x  / 0 
for  fix  = 0 
for  pu  / 0 
for  p„  =0 


**  = {|W  (I  + PyVPz2)  *}* 

When  the  velocity  of  the  fluid  relative  to  the  cable  is  two-dimensional,  in  the  xz  plane, 
Equations  (41a,  ft)  reduce  to 

KOI  <"*>=&  (43) 


kxl  ox  dZ 


which  gives 


<oJx>  Ax  + \Px\3l2!Px 


for  Px  ^ 0 
for  pz  = 0 


When  the  velocity  of  the  fluid  relative  to  the  cable  is  three-dimensional,  it  is  often  con- 
venient (see  Section  4.1)  to  assume  that,  for  the  purpose  of  obtaining  equivalent  mean  fluid 
velocity  components  <wz)  and  <«*v>,  the  x and  y components  of  the  relative  fluid  velocity  can 
be  treated  independently.  Then  Equation  (44)  may  be  used  to  obtain  <a>x>,  and  a corresponding 
equation  in  the  y direction  may  be  used  to  obtain  (<«„>.  Rewriting  Equation  (41a)  in  the  form 

»t 

IOx>|<7x(l  + (Oy)2/(<Tx)i)i  = kx|®x(l  + Oy2/oz2)*  dZ  (45) 

J 0 

it  can  be  seen  that  by  replacing  ox  and  au  by  (ax')  and  <a*>  respectively  only  in  the  square  root 
term  of  the  integrand,  this  equation  reduces  to  Equation  (43).  A similar  result  may  be  obtained 
from  Equation  (4lh)  for  the  y direction.  If,  over  the  range  of  integration,  the  variation  of 
(I  -f-  oy2/ox2)*  is  high,  the  approximation  may  be  poor.  It  should  be  noted  that  this  independent 
treatment  of  the  velocity  components  (referred  to  as  “independent  velocity  component  assump- 
tion") is  confined  to  the  calculation  of  equivalent  mean  fluid  velocity  components;  the  uniform 
velocity  three-dimensional  results  of  Section  3.3.2  should  be  used  to  obtain  angular  deviations 
in  pitch  and  roll. 

In  the  following  presentation  of  winds  above  the  sea  and  ocean  currents,  subscripts  A and 
W replace  subscript  F when  the  fluid  is  “air"  and  “water”  respectively. 


4.1  Winds  Above  the  Sea 

In  representing  the  wind  above  the  sea,  it  is  assumed  that  wind  direction  does  not  vary  with 
altitude  and  that  the  helicopter  heads  into  wind.  This  results  in  </ja  = <!>ao  = ±180°.  Generally, 


J. 
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conditions  above  the  sea  are  very  difficult  to  predict  and  a number  of  different  velocity  functions 
are  possible  for  various  stability  classifications.7  The  influence  of  the  rotor  downwash  compli- 
cates the  air  flow  even  further.  In  view  of  these  uncertainties,  a mathematically  complex  function 
is  not  warranted.  In  fact,  one  may  argue  that  a uniform  velocity  is  sufficient,  in  which  case  the 
results  already  obtained  in  Section  3.3  may  be  used  directly. 

Two  commonly  used  wind  profiles,  known  as  “power”  and  “log”  laws,  are  assumed  here 
and  are  defined  below.  The  log  law  profile  is  as  used  for  the  neutral  stability  classification.7  Each 
profile  requires  wz  0 at  the  sea  surface  and  wz  =-l  at  P0.  The  dimensionless  quantity  H is 
defined  as  h/zi,  where  h is  the  height  of  P0  above  the  sea  (i.e.  hover  height). 

(a)  Power  Law 

-wx=(\  -Z/H)"*  (46) 

where  n is  a positive  integer,  commonly  assumed  to  be  six. 

( b ) Log  Law 

- wz  — In  {c  — (e  — 1)  Z/H}/ln  c (47) 

where  c = 1 + h/ho  and  h0  is  the  surface  roughness  length  (on  p.  140  of  Ref.  7,  h0  = 
2-5  x 10  4 m for  calm  sea,  fto=  4 0 x 10  5 m for  rough  sea). 

The  above  two  formulae  give  very  similar  wind  profiles  provided  the  appropriate  values 
for  n and  c are  chosen. 

For  the  wind  profiles,  the  velocity  of  the  fluid  relative  to  the  cable  is  two-dimensional 
provided  any  helicopter  motion  is  in  the  same  direction  as  the  wind  (i.e.  fw=0or  ±180°, 
or  \y  0);  otherwise,  the  relative  velocity  is  three-dimensional.  Using  the  “equivalent  mean 
fluid  velocity”  concept,  for  the  two-dimensional  case,  (oiy)  = 0 and  <a>x>  is  readily  given  by 
Equation  (44).  However,  for  the  three-dimensional  case,  <eox>  and  <o» „>  should  be  recalculated 
for  each  value  of  \y.  It  is  therefore  convenient  to  use  the  independent  velocity  component 
assumption  outlined  in  Section  4.  Then,  for  a given  value  of  Ax,  the  values  (wx)  and  <tuy>  remain 
unchanged  for  all  values  of  \y.  For  the  wind  profiles  assumed,  using  typical  parameter  values 
experienced,  this  approximation  was  generally  found  to  be  good  when  |Ay|  was  reasonably  small 
(i.e.  less  than  about  0-5).  This  is  to  be  expected  since,  for  small  \y,  (1  ± oy2/ox2)*,  which  equals 
{1  + Aj,2/(coz  — Ax)2}4,  does  not  differ  significantly  from  unity  throughout  most  of  the  range  of 
integration  for  both  profiles. 

Analytical  solutions  for  <<ox>  from  Equation  (44)  are  possible.  However,  they  are  given  by 
lengthy  formulae  so  that,  except  for  special  cases,  results  presented  here  are  obtained  by  integrat- 
ing numerically  using  Simpson’s  formula.  For  the  special  case  when  A*  = 0 and  H = 1,  from 
Equation  (44),  Swx)  is  given  analytically  for  the  power  law  by 


■ 

* 1 

i 

■ 

(ojz)  — 

- i 

(1  — Z)2/n  dZ 

0 

= 

n 

2 + n 

and  for  the  log  law  by 

r .'•i 


— <«**>  = 


I In 2 {<■  — (c  - 

i In 2 c 

.J  o 


- I )Z} 


A Z 


* 

— 

2 

H 

c-  M 

H 

In 2 c 

(48) 


(49) 


Figure  1 1 shows  both  profiles  for  a range  of  parameter  values  n and  c.  Pairs  of  the  parameters 
have  been  suitably  chosen  to  give  similar  profiles  when  each  parameter  is  used  in  the  appropriate 
law.  At  A*  0 and  for  a number  of  values  of  \/H,  points  are  shown  on  the  curves  at  which 
<ai*>  = wz.  Using  the  same  parameter  values,  Figure  12  shows  the  ' ariation  of  <a>*>  with 
Ax  at  H I.  The  sudden  change  in  <<ux>  observed  for  each  curve  is  pi  Juced  by  the  effect  of  a 
reversal  in  the  fluid  drag  force  along  the  cable,  i.e.  <jz  changes  sign. 

Figure  13  shows  comparisons  between  two-dimensional  “general”  and  “approximate” 
solutions  for  0o  with  the  power  law  wind  profile  used.  The  former  solution  is  obtained  from 
Equation  (22a)  with  <f>  = 0 and  V = 0,  and  the  latter  is  given  by  Equation  (36)  with  the  equiv- 
alent mean  fluid  velocity  used.  As  with  the  uniform  velocity  solutions,  agreement  is  observed 
to  be  poor  where  large  angular  deviations  from  the  vertical  are  experienced.  Results  for  the  log 
law  were  found  to  be  very  similar. 
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4.2  Sea  Currents 

Like  winds  above  the  sea,  sea  currents  are  extremely  difficult  to  predict.8  One  may  well 
similarly  argue  that  any  function,  other  than  a uniform  one,  is  not  warranted,  in  which  case 
results  obtained  in  Section  3.3  may  again  be  used  directly.  However,  for  a non-uniform  velocity, 
the  classical  Ekman  spiral8  provides  a mathematically  convenient  function  that  is  fairly  realistic 
for  deep  water,  i.e.  particularly  for  oceans.  The  dimensionless  quantity  D is  defined  as  d/zi, 
where  d is  the  depth  of  frictional  resistance.  The  magnitude  and  direction  of  the  spiral  is  then 


defined  for  the  southern  hemispheref  and  with  4>a  = ±180°  by 

ut  = e ”Z,D  (50a) 

<fiw  = 3tt/4  - ttZ/D  (50 b) 

where  ipw  is  in  radians.  Equations  (48a,  b)  may  be  expressed  in  component  form  as 

wz  — in  cos  0H'  = e ”ZlD  cos  (3w/4  — irZ/  D)  (51a) 

oiy  = w sin  iJjw  = e ”Z,D  sin  (3jt/4  — irZ/D)  (516) 


The  above  function  is  illustrated  in  Figure  14.  At  the  sea  surface,  the  direction  of  the  current 
is  45"  to  the  left  (right,  in  the  northern  hemisphere)  of  the  wind  direction.  At  the  depth  of  friction, 
i.e.  when  Z D,  the  water  velocity  is  0 043  times  the  surface  velocity,  i.e.  o>  = 0 043.  The  sig- 
nificance attached  to  the  depth  d is  that  below  it,  the  velocity  is  considered  negligible. 

For  the  cable  length  above  the  sea,  the  velocity  uao  at  Po  is  generally  known  from  helicopter 
instrumentation.  For  the  submerged  cable  length,  the  water  velocity  uwo*  at  Po,  as  well  as  the 
depth  of  frictional  resistance,  d,  are  not  generally  given  by  any  instrumentation.  However, 
empirical  expressions  for  both  uw o*  and  r/have  been  derived.  From  References  8 and  9 respectively, 

d x 7-6  uao*  cosec*  </>e  (52a) 

uwo*  ~ 0-03  uao*  (52 b) 

where  </>e  is  the  geographical  latitude  and  d is  in  metres  for  uao*  in  metres  per  second.  At 
</>e  = 35",  Equation  (52a)  gives  d x 10  uao*.  The  above  references  assume  u^o*  is  measured  at 
an  altitude  of  15  m and  10  m respectively.  For  the  assumed  wind  profiles,  uao*  does  not  vary 
significantly  with  altitude  in  the  range  expected  for  a hovering  helicopter  with  a sonar  body 
suspended  by  a cable,  i.e.  roughly  10  m to  30  m.  Hence,  adjustment  to  the  above  formulae  is 
not  warranted  for  changes  in  altitude  within  this  range. 

In  Reference  8,  at  wind  velocities  below  about  6 m/s,  Equation  (52 a)  is  replaced  by 

d x 3-67 (u ao*)3  2 coscc*  </>e  (52 c) 

which,  at  4>e  35",  gives  d x 4-8(w*o*)3'2.  However,  this  further  refinement  is  not  considered 

necessary  in  view  of  the  limited  confidence  in  the  model  of  the  sea  current.  Because  the  angular 
deviation  of  the  cable  is  roughly  proportional  to  the  square  of  the  velocity  (see  Equation  (36)), 
errors  are  minimized  over  a large  range  of  velocity  by  choosing  the  empirical  relation  given  by 
Equation  (52a),  which  is  more  appropriate  for  higher  velocities.  Neither  relation  allows  a mean- 
ingful estimate  for  d at  the  equator,  i.e.  at  <\>R  — 0. 

During  normal  sonar  operation,  it  is  generally  desirable  for  the  sonar  body  to  reman 
vertical.  For  the  mathematical  model  used  here,  this  requires  (for  yi  = 0)  the  relative  fluid 
velocity  at  the  body  to  be  zero,  i.e.  A*  wz\  and  A„  = wyi,  where  wzi  and  cu„i  are  the  values  of 
wz  and  uiv  respectively  at  P\.  Equivalent  mean  fluid  velocities,  which  are  obtained  by  integrating 
numerically  using  Simpson’s  formula,  are  therefore  presented  assuming  this  condition. 

In  Figure  15,  both  the  Ekman  spiral  and  equivalent  mean  fluid  velocity  are  shown  in  the 
xy  plane  as  functions  of  Z/D  and  I/D  respectively.  Equivalent  mean  fluid  velocity  values  cal- 
culated using  the  independent  velocity  component  assumption  are  compared  with  the  values 
obtained  without  this  assumption,  i.e.  using  Equations  (42a,  b).  Agreement  is  seen  to  be  very 
good  for  the  lower  and  upper  range  of  values  of  I/D  shown.  As  explained  in  Section  4,  regions  of 
good  agreement  imply  there  is  only  small  variation  in  the  term  (1  + o'*2/®*2)*  throughout  the 
range  of  integration.  The  absolute  error  in  the  approximate  values  is  given  by  the  distance 

t For  the  northern  hemisphere,  </>w  = — 3»r/4  + nZ/D. 
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between  points  on  the  two  mean  velocity  curves  of  Figure  15  which  have  the  same  value  of 
1 1 ID.  This  error  is  a maximum  of  0 037  at  \/D  » 1 05  and  is  about  0 01  at  1 /D  = 0-6  and  14. 

| Although  these  errors  are  much  larger  than  those  generally  experienced  when  considering  winds 

, above  the  sea,  they  may  well  be  acceptable  when  viewed  in  relation  to  uncertainties  already 

, existing. 

A three-dimensional  “approximate”  solution  for  60  and  <j> o is  obtained  from  Equations 
(34 a,  b)  on  using  the  equivalent  mean  velocity  values  shown  in  Figure  15.  For  appropriate  values 
of  the  parameters  p and  (/wo*,  solutions  of  angular  orientation  in  pitch  and  roll  are  given  in 
Figure  16  and  are  compared  with  solutions  obtained  using  the  “independent  component” 
equivalent  mean  fluid  velocities.  Errors  are  seen  to  correspond  to  those  observed  in  Figure  15. 
Results  obtained  from  the  “general”  solution  are  not  shown  since  they  were  found  to  agree 
i very  closely  with  the  “approximate”  solution  results  (i.e.  full  lines  of  Figure  16).  Good  agreement 

’ is  to  be  expected  because  the  angular  deviations  from  the  vertical  are  small. 

? 4.3  Example  of  Sonar  Body  Suspended  from  a Helicopter 

So  far,  solutions  have  been  presented  only  for  a length  of  cable  totally  immersed  in  a single 
Newtonian  fluid.  However,  the  main  case  of  interest  in  this  study  is  of  a sonar  body  suspended 
from  a helicopter.  Here,  the  cable  passes  through  two  fluid  layers  (air  and  water)  and  the  method 
of  approach  outlined  in  Section  2.1  should  be  used.t  Using  a power  law  for  the  wind  above  the 
sea,  and  an  Ekman  spiral  for  the  sea  current,  an  example  is  given  below  in  which  cable  orien- 
tation and  shape  are  determined. 

Constants  that  are  independent  of  the  fluid  are  g = 9-81  m/s2,  dc  = 20  mm,  Cd  — 1-2, 
and  uh*  = 0. 

The  assumed  conditions  in  air  are  pa  = 1 -226  kg/m3,  mg  = 6-2  N/m,  zi  = h = 20  m 
(i.e.  H = h/zi  = 1),  uao*  = 20  m/s,  <pA  = ±180°,  and  n = 6 (for  power  law). 

The  assumed  conditions  in  water  are  pw  = 1023  kg/m3,  T\  = 60  N,  z\  = 40  m, 
d = 10  uao*  = 200  m (i.e.  D = d/zi  — 5),  uwo*  = 0 03  uAo*  = 0-6  m/s,  y\  = 20°,  *i  = i/n  = 
(0-75jt  — 7 t/D)  rad  = 99°  (i.e.  at  attached  body,  cable  is  inclined  in  a vertical  plane  in  the  flow 
direction),  and  mg  = 30  N/m  (obtained  by  subtracting  the  upthrust,  pw  g (Jm/C)2,  from  the 
value  of  mg  in  air). 

In  water,  using  Equation  (21),  7o  = 7)  ± mgz\  = 180  N,  which  is  equal  to  T\  in  air.  The 
dimensionless  quantities  p and  U o*  are  then  calculated  for  both  air  and  water  and  are  given  in 
Table  2 together  with  a summary  of  the  parameter  values  required  in  the  calculations. 

On  numerically  solving  the  general  equations,  cable  orientation  and  linear  displacement 
coordinates  are  obtained.  Table  3 presents  these  values  at  the  sonar  body,  sea  surface,  and 
helicopter.  All  displacement  coordinates  are  expressed  in  terms  of  dimensionless  coordinates 
for  the  cable  in  air  and  should  be  multiplied  by  20  to  give  displacement  in  metres.  Figure  17 
illustrates  the  cable  shape  together  with  the  horizontal  fluid  velocity  vectors.  For  each  fluid,  these 
vectors  represent  the  scaled  velocity  w,  where  a»o  = 1 . 


t The  method  could  easily  be  extended  to  any  number  of  fluid  layers. 
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TABLE  2 

Parameter  Values  that  are  Dependent  on  Fluid 


Quantity 

Value  in 

water 

Value  in 
air 

PF 

1023  kg/m3 

1 -226  kg/m3 

r, 

60  N 

180  N 

mg 

3-0  N/m 

6-2  N/m 

Zl 

40  m 

20  m 

UFO* 

0 • 6 m/s 

20  m/s 

kc  = kpF  dr  Cd 

12-3  kg/m2 

0 0147  kg/m2 

(kr  Zl /7i)* 

2 • 86  s/m 

0-0404  s/m 

p = mgzi/Ti 

20 

0-69 

U0*  = Ufo*  = (kc  Zi/7-,)*  ufo* 

1 -72 

0-81 

TABLE  3 

Orientation  Angles  and  Displacement  Coordinates 


Quantity 

(1)  Sonar  body 

Value  at: 

(2)  Sea  surface 

(3)  Helicopter 

6 

-3-3° 

-17-5° 

-26-3° 

* 

-19-7° 

-29-0° 

-14-5° 

X 

-0-681 

-0-393 

0 

Y 

1-311 

0-404 

0 

Z 

3 

1 

0 

5.  CONCLUDING  REMARKS 

A three-dimensional  steady  state  mathematical  model  of  a cable  used  to  suspend  a sonar 
body  from  a helicopter  has  been  derived.  The  effect  of  the  sea  waves,  which  would  make  the 
problem  unsteady,  is  neglected.  Results  are  presented  in  dimensionless  form  for  a length  of  cable 
totally  immersed  in  a Newtonian  fluid,  and  may  therefore  be  applied  to  a number  of  practical 
problems. 

In  order  to  obtain  analytical  solutions  giving  cable  orientation  at  its  suspension  point,  it  is 
assumed  that  the  orientation  angles  are  small  and  that  the  cable  mass  is  concentrated  at  the 
attached  body.  Solutions  of  the  general  and  approximate  equations  are  obtained  for  the  case 
when  the  fluid  velocity  is  uniform.  In  general,  for  both  two-dimensional  and  three-dimensional 
cases,  where  maximum  angular  deviations  from  the  vertical  are  small,  agreement  between  the 
“general”  and  “approximate”  solutions  is  good ; large  deviations  usually  result  in  poor  agreement. 

The  use  of  an  equivalent  mean  fluid  velocity  to  replace  any  non-uniform  profile  allows 
solutions  to  be  readily  obtained  from  the  uniform  velocity  results  presented  in  the  figures.  This 
enables  a more  compact  presentation  of  results  for  the  various  non-uniform  fluid  velocity  pro- 
files used  here.  These  non-uniform  profiles  include  two-dimensional  “power"  and  “log”  laws 
for  the  wind  above  the  sea,  and  a three-dimensional  Ekman  spiral  for  the  sea  current.  The  above 
concept  would  also  be  very  useful  for  simplifying  the  representation  of  the  environment  in 
developing  a mathematical  model  for  a flight  simulator.  Comparisons  between  non-uniform 
fluid  velocity  solutions  of  the  general  and  approximate  equations  result  in  similar  conclusions 
to  those  found  for  the  uniform  fluid  velocity  case. 
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NOTATION 


Cd  Drag  coefficient  of  cable  in  the  lateral  direction  of  flow. 

D d/zi 

{Fz,  Fy,  Ft)  Fluid  force  vector 

[Gz,  Gy,  Gt)  Gravitational  force  vector 

H h/zi 

P Upper  end  of  cable  element 

Po  Upper  end  of  cable  length 

Pi  Lower  end  of  cable  length 

S s/zi 

T Tension  force  at  P 

{ Tz , Ty,  Tt)  Tension  force  vector 

U*  Non-dimensional  velocity  magnitude  (kc  zi/Ti)*  u* 

U,  V x,y  components  of  U* 

X,  Y,  Z x/zi,  y/zi,  z/zi 


a 

b 

c 

d 

dc 

g 

h 

h0 

kc 

m 

n 

<7 

s 

t 

u* 

UF*,  UH* 
U,  V,  W 
UF,  VF 


pIWoWoU 

1 + Ml  - Z) 

1 + h/ho 

Depth  of  frictional  resistance 
Cable  diameter 

Gravitational  acceleration  (along  zg  axis) 
Height  of  Po  above  the  sea  (i.e.  hover  height) 
Surface  roughness  length 
ipF  dc  Cd 


Effective  cable  mass  per  unit  length  (allowing  for  buoyancy  effect) 
Positive  integer 
(1  + 4 /a*)* 


Distance  of  P from  Po  along  cable 
sin  8 

Velocity  iff  magnitude  fluid  relative  to  cable 
Velocity  magnitude  of  fluid,  cable  relative  to  earth 
x,  y,  z components  of  velocity  of  fluid  relative  to  cable 
x,  y components  of  fluid  velocity  in  earth  axes 


i 


j 


X 


UH,  VH 

xyz 
x'y'z' 
xeyezn 
x,y,z, 
x,  y,  z 
0,  <D 
0/,  O, 

«.  P,  y 

*x 

“v 

Px 

P« 

8 

o,  <f> 

A 

Ax,  Av 

P 

PF 

a 

Qx%  Gy 
<t>E 

X.  Y 
<!> 

<I>F,  <pH 
to 

WZ,  Wy 
{} 

O 

Subscripts: 
A,  W 
s 

0,  1 


x,  y components  of  cable  velocity  in  earth  axes 
Cable  axes 

I^'ermediate  rotation  axes 

Earth  axes 

Cable  element  axes 

Linear  displacement  components  of  P in  cable  axes 
Matrices  representing  rotations  through  0,  <f> 

Matrices  representing  rotations  through  0 + 80,  <f>  + &</> 

Direction  cosine  angles  of  cable  element  axes  with  respect  to  cable  axes 

mo  + ISvW)-*}* 
mo  +/3xW)-*}+ 

[ U*UdZ 


J o 


Denotes  increment  in  quantity  prefixed 

Pitch,  roll  Euler  angles  when  rotating  xyz  axes  to  x,y,z,  axes  (pitch  followed  by 
roll) 

Uh*/Ufo* 

x,  y components  of  A 

mgzi/Ti 

Fluid  density 

U*/Ufo* 

x,  y components  of  a 
Geographical  latitude 

Yaw,  pitch  Euler  angles  when  rotating  xyz  axes  to  x,y,z,  axes  (yaw  followed  by 
pitch) 

Velocity  direction  of  fluid  relative  to  cable 
Velocity  direction  of  fluid,  cable 
Uf*/Ufo* 

x,  v components  of  w 

Used  to  write  a column  vector  as  a row  vector 
Enclosed  quantity  is  a mean  value 


Replaces  subscript  F when  the  fluid  is  air,  water 

Quantity  is  in  cable  element  axes 

Value  at  Po,  Pi  (may  follow  other  subscripts) 
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APPENDIX  A 

Cable  Shape  from  General  and  Approximate  Equations 

(1)  General  Equations 

In  dimensionless  form,  derivatives  of  the  cable  coordinates  X,  Y and  cable  length  S with 
respect  to  Z are  given  from  Equations  (la,  b,  c)  by 

dX/dZ  = tan  8 

d Y/dZ  = — sec  9 tan  <f> 

dS/dZ  = sec  8 sec  <f> 


On  integrating  between  P and  Pi, 


Xi  — X = 


tan  8 dZ 


• i 


Yi  — Y = — J sec  8 tan  <f>  dZ 


I 


Si  — S = | sec  8 sec  <f>  dZ 
2 


The  left-hand  side  of  each  of  the  above  equations  is  first  determined  at  all  required  equispaced 
values  of  Z by  integrating  numerically  for  Z decreasing  from  one  to  zero.  At  Z = 0,  the  variables 
X,  Y,  and  S aie  zero  so  that  Xi,  Yi  and  Si  are  given  directly.  This  enables  X,  Y,  and  S to  be 
obtained  for  0 < Z < 1 from  the  previouslyscalculated  left-hand  sides.  The  above  integrals  are 
calculated  using  Simpson’s  formula  following  each  step  of  the  numerical  method  used  to  obtain 
8 and  <f>.  This  procedure  is  particularly  suited  to  the  Runge-Kutta  second  order  method  used 
(see  Section  3.1.1),  since  the  mid-interval  values  required  by  Simpson’s  formula  are  directly 
provided  by  the  Runge-Kutta  method. 


(2)  Approximate  Equations 

For  small  angles,  the  above  equations  may  be  written  ast 


Xi  — X x 


8 dZ 


Yi 


<f,dZ 


/•< 


Si  — S x 


[1  + W2  + <f>2))  dZ 


t Second  order  terms  in  8 and  4>  are  retained  for  the  expression  for  Si  — S',  otherwise 

SxZ. 


Generally,  for  a variable  fluid  velocity,  the  procedure  for  obtaining  X,  Y,  and  S is  as  described 
above  in  Part  (I).  For  a uniform  fluid  velocity,  or  a non-uniform  one  if  the  equivalent  mean 
fluid  velocity  components  are  obtained  analytically,  0 and  <f>  are  given  analytically  as  linear 
functions  of  Z (see  Equations  (34a,  b )).  On  substituting  these  functions,  integrating  analytically, 
and  then  using  Equations  (34a,  b)  to  simplify  the  result,  we  have 

Xi-  X=i(l  - Z)(0  i + 0) 

Yi  - y = -i(l  - Z)(<h  + <f>) 

Sx  - 5 = (1  - Z)[l  + m2  + »i9  + 02  + *i*  + M + **)] 

The  variables  X,  Y,  and  S can  be  obtained  as  before,  but  because  6o  and  4>o  are  given  directly 
without  any  numeiical  integration,  use  can  be  made  of  these  values  at  commencement  of  the 
calculations.  At  Po,  the  variables  X,  Y,  Z,  and  5 are  all  zero;  hence  the  values  X\,  Yi,  and  Si  are 
given  by 

Xi  = |(0i  + 0o) 

Yi  = — i(^i  + 4>o) 

Si  = 1 + i (0i2  + 0i0o  + 0o2  + ^i2  + ^i^o  + ^o2) 

Using  the  two  above  sets  of  equations  together  with  Equations  (34a,  b),  the  variables  X,  Y,  and 
S are  given  directly  for  any  value  of  Z. 
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Fig. 3 Comparison  between  two-dimensional  'general'  (full 
line)  and  'approximate'  (broken  line)  solution  of 
angular  orientation  for  uniform  fluid  velocity; 

60  vs.  UQ  for  parameters  \i  and  0i 
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Fig. 4 As  Fig. 3,  but  0O  vs.  p for  parameters  6i  and  U( 
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Fig. 6 Comparison  between  cable  shapes  obtained  using 
'general'  (full  line)  and  'approximate'  (broken 
line)  two-dimensional  solution  for  uniform  fluid 
velocity 


a)  independent  variable  i|»Q  (deg) 
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b) independent  variable  X\ (deg) 
Yi  = 10° 
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Fig. 7 Comparison  between  three-dimensional  'general' 

(full  line)  and  'approximate'  (broken  line)  solution 
of  angular  orientation  in  pitch  and  roll  for  uniform 
fluid  velocity,  parameter  y,  and  independent  variable 
(a)^0  and  (b)xi«  In  Figs. 7 to  9,  tick  marks  indicate 
equispaced  values  of  the  independent  variable;  values 
are  shown  only  in  Fig. 7a,  but  are  the  same  at 
corresponding  circumference  positions  for  each  figure 
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Fig. 10  Comparison  between  cable  shapes  obtained  using 
'general'  (full  line)  and  'approximate'  (broken 
line)  three-dimensional  solution  for  uniform  fluid 
velocity.  Four  views  of  same  shapes  are  shown 
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Fig. 11  Velocity  profiles  for  wind  above  the  sea;  power  law 
(full  line)  for  parameter  n,  and  log  law  (broken 
line)  for  parameter  c.  Points  at  which  <ojx>  = wx 
are  shown  for  Xx  = 0 and  parameter  1/H 
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Fig. 13  Comparison  between  two  dimensional  'general'  (full 
line)  and  'approximate'  (broken  line)  solution  of 
angular  orientation  for  power  law  with  equivalent 
mean  fluid  velocity  used  in  latter;  0O  vs.  Xx  for 
parameters  y and  Oj 
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Fig. 14  Ekman  spiral  velocity  profile  for  sea  current 
(viewed  below  sea  level) . Horizontal  velocity 
vectors  are  shown  at  intervals  in  Z/D  of  0.1 
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Fig.  15  Ekman  spiral  (u)x,u)y)  and  equivalent  mean  fluid 
velocity  (<iox>,<u)y>)  for  independent  variables 
Z/D  and  1/D  respectively  (values  shown  beside 
tick  marks) . Equivalent  mean  fluid  velocity 
calculated  using  'independent  velocity  component 
assumption'  is  also  shown 
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Fig. 16  Three-dimensional  'approximate'  solution  (full 
line)  of  angular  orientation  in  pitch  and  roll 
for  Ekman  spiral  with  parameter  y and  Uwo*  and 
independent  variable  1/D  (values  shown  beside 
tick  marks).  Solution  is  also  shown  (broken 
line)  using  'independent  velocity  component 
assumption ' 
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Fig. 17  Example  of  sonar  body  suspended  from  a helicopter; 

note  the  difference  in  scales  of  the  fluid  velocity 
vectors  for  air  and  water 
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